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ABSTRACT 

We present a description of the CL-based package xdspres, which aims at being a 
complete reducing facility for cross-dispersed spectra taken with the Ohio State Infrared 
Imager/Spectrometer, as installed at the SOAR telescope. This instrument provides 
spectra in the range between 1.2|j,m and 2.35iJ.m in a single exposure, with resolving 
power of R ~ 1200. xdspres consists of two tasks, namely xdflat and doosiris. The 
former is a completely automated code for preparing normalized flat field images from 
raw flat field exposures. Doosiris was designed to be a complete reduction pipeline, 
requiring a minimum of user interaction. General steps towards a fully reduced spectrum 
are explained, as well as the approach adopted by our code. The software is available 
to the community through the web site http://www.if.ufrgs.br/^ruschel/software. 

1. Introduction 



X 



Cross-dispersed spectroscopy makes possible to acquire information of wide spectral regions 
in a single exposure, by projecting several dispersion axes on the detector simultaneously. As 
a consequence, the reduction process required to analyze this kind of data is complicated, since 
different diffraction orders need to be selected, extracted, calibrated independently and combined 
in the final step. This difficulty led many authors to dev elop methods and software packages for 



the reduction of cross-dispersed and echelle spect ra (e.g. iMoreno et al. Ill982l : iRossi et al. Ill985l : 



Piskunov &: Valenti 



2002 



Bochanski et al. II2009I ). 



In the past decade t he near infrared (NI R) has also been explored by cross-dispersed spec- 



trographs, such as Spex (JRavner et al. 1 120031 ) at the NASA Infrared Telescope Facility (IRTF), 
with a resolving pow'er of '^ 2000 and reaching from 0.8 to 5.5|J.m. Other examples are Triple - 
Spec (JEdelstein et ahlboOTl l and the Folded-port Infrared Echellette (FIRE) (jsimcoe et ahlbood ). 
achieving R ~ 2600 and R ~ 6000 respectively, and covering roughly the same wavelength domain 
(0.8 - 2.4|j.m). 



Another instrument of similar capabilities is the Ohio State Infrared Imager/Spectrometer 
(OSIRIS), currently installed at the Southern Astrophysics Research Observatory (SOAR), attached 
to the 4.1m telescope. OSIRIS provides spectral coverage from 1.0|J.m to 2.4|j,m in cross-dispersed 
mode, with a resolving power of ~ 1200. High resolution (R ~ 3000) long-slit modes are also 
available, but multi-band spectroscopy of this kind suffers from differences in aperture and seeing. 

However, reduction of NIR spectra has a complexity of its own, mostly related to telluric 
spectral features, both in absorption and emission, and black body radiation due to the telescope 



itself . A rich literature has been developed on the subject (e.g. iMaiolino et al. lll996l : IVacca et al 



20031 : ICushing et al. Il2004l l. 



There are currently no specific software packages available for the reduction of cross-dispersed 
spectra taken with OSIRIS. Aiming at providing a fast and highly automated task, we developed 
the XDSPRES (acronym for cross-dispersed spectra reduction script) package. The CL language was 
chosen due to the availability of almost all of the basic tasks needed to perform the reduction in 
the Image Reduction and Analysis Facility (IRAF) software (JTodv Ill986l . Il993l ). 



In ^ we describe main aspects of the instrument, focusing on its effects on the reduction 
process. In ^ we describe the general steps towards a fully reduced spectrum, as well as the 
approach adopted by the xdspres package to each of these steps, and finally in §1] we give a brief 
summary. 



2. OSIRIS 

In this section we discuss the main aspects of the cross-dispersed mode of OSIRIS, with special 
attention to those characteristics that are relevant to the reduction process. A complete description 
of the instrument can be found in its on line User's Manuao 



The detector is a 1024x1024 HAWAII array (JHodapp et al. Ill996l ). sensitive to wavelengths of 
up to 2.5|xm. Equation [1] models the non-linear behavior of the array, which only becomes critical 
above 28,000 counts. Usually the detector is read only at the end of the integration, but since it 
can be read non-destructively different sampling methods could be implemented. 

A residual image is sometimes seen, specially when bright sources are observed in acquisition 
mode. This means that eventually some of the first spectra taken after the target acquisition images 
have to be discarded. Residuals have approximately 2% of the intensity of the original source, and 
it should not be a problem to science exposures that have typical counts below one thousand. 



^http://www.ctio. noao.edu/instrunients/ir .instruments/ 
osiris2soar/manual/ 



ADU' 

— — - = 1.00108 - 1.015777 x lO^^ADU 

+ 1.548099 X lO^^^ADU^ 

- 1.945376 X lO^^^ADU^ (1) 



In cross-dispersed mode OSIRIS projects almost six orders on the detector, from which three 
are extracted. Wavelength coverage for each of the extracted orders are 1.2 - 1.5, 1.5 - 1.9 and 1.9 
- 2.35|J.m, for the J, H and K bands respectively, all of them with R ~ 1200. Orders that are not 
extracted include a small portion of the J band (1.0 - 1.2|j,m), and second order duplicates of the 
J band, located to the right of the K band. Figure [1] shows an example of sky spectrum where the 
three main orders are evident. Orders that are not extracted are also visible in figure [2j 

From figured] it can also be noted that dispersion axes are nearly vertical, meaning that within 
a given aperture each line corresponds to a particular wavelength. The misalignment between 
detector lines and wavelength coordinate are less than one pixel from one end of the slit to the 
other, or less than one third of the full width at half maximum (FWHM) of a emission line in the J 
band. Therefore corrections to dispersion axis orientation were not attempted, and all extractions 
assume a vertical dispersion. 



3. Reduction process 

3.1. Flat field 

In cross-dispersed mode, flat fleld images are taken with the cross-dispersing grism already 
positioned, which results in a spectrum of the flat field lamp, rather than an evenly illuminated 
image. Since the main purpose of a flat field is to identify pixel-to-pixel variations which are 
intrinsic to the detector, the continuum that corresponds to the spectral energy distribution of the 
lamp has to be removed. Moreover, two sets of flat fields are needed, one with the fiat field lamp on 
and another with the lamp off. The later is required because thermal radiation from the telescope 
becomes appreciable in the low energy end of the spectrum, as it can be seen in Fig. [2j Typical 
sets consist of 10 exposures of each kind. 

The xdflat task automates the preparation of a normalized fiat field image, which will be later 
used to correct the science images. First it applies a linearity correction to all flat field images, 
according to equation [H Then both sets (fiat-on and flat-off) are averaged independently, and 
the resulting flat-off image is subtracted from the flat-on. We have omitted a figure showing the 
subtracted flat because it is visually identical to the flat-on. The only noticeable difference is the 
suppression of a few hot pixels at the lower portion of the image. 

To remove the spectrum of the flat field lamp xdflat begins by extracting each order. Apertures 



are identified by a centering algorithm (apfind) that searches for three local maxima in the central 
lines of the chip. The peaks are assumed to be separated by more than 30 pixels and have an 
approximate width of 80 pixels. Aperture sizes are reevaluated by setting the borders at 20% 
of the peak intensity of each order. A tracing algorithm (aptrace) moves in regular five pixel 
steps along the dispersion axis, assessing changes in peak location for each order, leading to a two- 
dimensional description of the aperture position. The aperture tracing function, a second order 
Legendre polynomial, is fitted to predefined sample regions of the chip that are less affected by 
scattered light. Errors in the two-dimensional aperture border definitions are usually below three 
pixels. 

A 30th order Legendre polynomial is fit to the spectrum, which is then normalized. Such a 
high order polynomial is justified by the complex pattern produced by the fiat field lamp as it passes 
through the spectrometer, as shown by figure El Artificial oscillations at the apertures' limits are 
ignored after extraction. Typical RMS of the fit is below 5000 ADU, which may seem high but 
actually amounts to roughly 2% of the average signal. The final fiat-field image has all its pixel 
counts set to 1, except those on the regions occupied by the spectrum, which are replaced by the 
ratio between the original count and the fitted polynomial. 



3.2. Subtraction Object - Sky 

In the NIR spectral region the atmosphere plays an important role. Besides a significant telluric 
absorption, several atmospheric emission lines are entangled with the spectrum of the astronomical 
source (see figured] for a sample spectrum of the sky, where the J, H and K bands are identified). 

The process of removing telluric emission lines is commonly known as sky subtraction, or sky 
chopping, and the angular size of the target dictates whether additional off-source exposures are 
required. In the case of point sources, which occupy only a small fraction of the slit, one can take 
exposures with the source in two different positions along the slit, and later subtract subsequent 
images. This is the technique employed to obtain the spectra of standard stars, and it makes more 
efficient use of telescope time. When extended sources are concerned, a separate set of exposures 
taken from a nearby dark region of the sky is needed, a process commonly referred to as nodding. 

The doosiris task was developed to reduce spectra from extended sources, therefore it assumes 
that a set of sky exposures was taken along with the science exposures, in order to remove the 
telluric emission lines. There are two ways by which users can inform the software about the 
nature of each image, namely: interactively identifying them via SAO Image DS9, or providing 
an ASCII file with the type of exposure with respect to its numerical order. For further details 
refer to the xdspres Manual. No attempt was made to provide a software solution for identifying 
different types of exposures, as specific criteria regarding the spectrum of the astronomical target 
would have to be predefined, adding, in our judgement, unnecessary complexity to the code. 

Nodding patterns that make best use of telescope time use each sky exposure in more than one 



subtraction, as in 0-S-O or 0-S-O-O-qj. It is thus impractical to simply subtract a combination of 
sky images from an equivalent combination of target ones. Instead of assessing the relevant physical 
quantities, a routine searches for the best telluric calibrator based on the file name index, assuming 
that these are sequentially numbered after the time of exposure. 



3.3. Extraction and Wavelength Calibration 

Extraction of science spectra follows the same procedures that were described in section 13. Jn . 
The sky spectrum is extracted using the same aperture definitions of the target spectrum. 

Wavelength calibration is based on strong OH lines present in the sky exposures, a sample of 
which is shown in figure [5j As of the moment of the publication of this paper, OSIRIS presents 
what appears to be an illumination problem that produces lines across the detector, in the direction 
perpendicular to the dispersion axis. Since these lines can lead to confusion in the OH line iden- 
tification process, a high order polynomial is used to fit and remove the vertical profile identified 
between columns 980 and 1024 (see figure H]). 

Interactive line identification is usually the best option, and since the dispersion function is 
almost linear there is no need for manually identifying more than four well spaced features. If the 
dispersion function fitting was successful, the unid entified features will ma tch those in the line list 



provided with xdspres, which was extracted from lOliva fc Origlia I (Il992l ). Doosiris also provides 
an option to automatically identify OH features in the spectrum of the sky that uses the reidentify 
task, which requires a previously identified image. For ~ 20 identified features, typical residuals 
are below 2 pixels, which translates into roughly it 50 km s~^. 



3.4. Telluric Removal and Flux calibration 

The subtraction of sky exposures from the on source images, obviously can only account for 
emission features of the atmosphere. To deal with the more subtle problem of removing atmospheric 
absorption doosiris by default uses the spectrum of an AOV star, that should be obtained just before 
or after the science images. If the observed standard star has a different spectral type, the model 
atmosphere spectra, mentioned below, have to be replaced accordingly. 

The standard star, being a point source, do not need a separate set of sky exposures, because 
it occupies only a small fraction of the slit. Doosiris is prepared to manage two or three different 
star positions on the slit. In either case subsequent exposures are subtracted and the resulting 



^Where "O" stands for object and "S" for sky 

'Although doosiris is prepared to automa 
spatial profiles should be personally reviewed 



'Although doosiris is prepared to automatically define the borders of each aperture, targets that have complex 



images are summed; a sample of this sum can be seen in figure [6l After division by the normalized 
flat field image, both spectra are extracted and summed. 

The spectrum of an AOV is almost devoid of metallic absorption lines, but the H lines that 
are present need to be eliminated before it can be applied to t he science spectruni as a telluric 



calibrator. The method employed here follows the reasoning of IVacca et al. I (|2003l ). but with a 
different implementation. It basically consists in dividing the spectrum of the standard star by a 
model atmosphere of Vega, obtained from R. Kuruc^. 

First the model of Vega was smoothed by a Gaussian with a equal to the FWHM measured 
in a NeAr calibration lamp, to match the resolving power of the standard star. A spline was then 
adjusted to the continuum, leading to a purely absorption spectrum. The later is provided with 
the XDSPRES package. The actual division of the reference star is performed by the telluric task, 
which allows for the shifting and scaling of the model. Figure [7] shows a comparison between the 
observed spectrum of the standard star and a model atmosphere for Vega. 

Once the absorption lines due to the stellar atmosphere have been removed, the spectrum 
becomes essentially a black body with telluric features. Its normalization by a polynomial that 
acts as a pseudo-continuum returns a purely telluric spectrum. Unabsorbed regions, that translate 
into sample regions for continuum fitting, were identified with the aid of NSO/Kitt Peak FTS data 
produced by NSF/NOAQj. This division of the science spectrum also allows shifting and scaling. 
Some of the strongest telluric bands cannot be fully removed. Additionally, the high absorption in 
these regions causes a significant decrease in S/N. 

The same polynomial employed as a pseudo-continuum for the reference star is later used to 
produce an independent sensitivity function for each aperture, by comparing it to a black body 
of 9480 K. This procedure restores the correct slope of the spectrum regardless of the accuracy in 
absolute flux. The later is estimated from the exposure time and magnitude of the standard star, 
which has to be provided by the user. Figure [8] shows the effects of telluric line removal and flux 
calibration to a sample spectrum. 

One would expect that a good flux calibration leads to a perfect alignment of the spectrum 
between different apertures. Although generally true, it has been observed that agreement is harder 
to achieve where the H and K bands meet. Strong telluric absorption bands near 1.9|xm difficult 
the evaluation of the sensibility function causing large deviations in the final spectrum. Figure [9] 
shows a completely reduced spectrum encompassing the whole spectral range. 



''http://kurucz.harvard.edu/stars.html 

^Available at http://www.eso.org/sci/facilities/par anaI7] 
instruments/isaac/tools/spectroscopic_standards.html 



4. Summary 

We have presented the xdspres CL-based package, consisting of the xdflat and doosiris tasks, 
aimed at being a complete reduction facihty for cross-dispersed spectra taken with the OSIRIS 
spectrometer, currently installed at the SOAR telescope. This particular instrument provides a 
relatively large spectral coverage, being able to project the full range between 1.2|a.m and 2.35|xm 
over the detector in a single exposure. The blazing of different orders in the same image adds 
complexity to the already lengthy reduction of infrared spectroscopy data, xdspres automatically 
performs the more mechanical and time consuming steps of the reduction, at the same time that 
it allows considerable user interaction in the more subjective stages. In addition, the possibility of 
a fast reduction provides means to make site adjustments to the observation s trategy. As a sampl e 



of actually published data that was fully reduced with the xdspres tasks, see lRiffel et al. I (|201ll ). 
The complete software package and its documentation is available to the community at the web 
site http://www. if. ufrgs. hr/^ruschel/ software. 
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Fig. 1. — Spectrum showing atmospheric emission hnes and identifying the orders that are projected 
on the chip. The exposure time for this image was 20s. Horizontal hnes seen in this and subsequent 
images, mainly around 140 and 650 pixels, are probably the result of scattered light that reached 
the grating. 
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Fig. 2. — Left A sample flat field image with the lamp turned ofl'. Thermal radiation from the 
telescope can be seen, as the flat field exposure is taken with grism already positioned. Right A 
flat fleld with the lamp turned on, clearly showing the three orders in the center. Also visible 
are: further J band orders at the lower left and beyond the K band to the right; hot pixels in the 
lower corners; two groups of cold pixels near the center of chip; a small portion of an order at the 
detector's left border, between lines 400 and 800. Both images were taken with 3.2s of exposure. 
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Fig. 3. — Mean flat fleld spectrum at the J band (grey), and fitted function (black), RMS for this 
fit is ~ 4600 ADU, which correspond to roughly 2% of the average signal. Artificial oscillations 
of the fit have no practical effects over the science spectrum as these portions are ignored after 
extraction. 
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Fig. 4. — Average of chip columns 980 to 1024. This background profile is removed prior to 
wavelength calibration to avoid confusion with OH emission lines. 
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Fig. 5. — A sample of the sky's spectrum at the H band. 
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Fig. 6. — The resultant image of subtracting subsequent exposures of the standard star occupying 
two different positions along the slit. 
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Fig. 7. — Comparison between the standard star H band spectrum (top) and a model atmosphere 
(bottom) for Vega. The resolving power of both spectra is 1200. 
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Fig. 8. — From top to bottom: a) Sample spectrum at J band just after wavelength calibration, 
b) Same spectrum after the removal of telluric lines, c) Flux calibrated spectrum. Areas of strong 
telluric absorption cannot be fully corrected, and even if they could the signal to noise ratio would 
still be much lower than the rest of the spectrum. 
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Fig. 9. — Flux calibrated spectrum in the whole spectral range. Aperture transitions are at 1.5|am 
and 1.9|J.m. Strong telluric absorption bands that dominate the spectrum between 1.8 and 2.0|a.m 
difficult the alignment between the H and K bands. 



